Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





The proper motion of the isolated neutron star 

RXJ1605.3+3249* 

C. Motch\ K. Sekiguchi^, F. Haberl^, V.E. Zavlin^'^, A. Schwope^, and M.W. Pakull^ 

^ Observatoire Astronomique, UMR 7550 CNRS, 11 rue de I'llniversite, F-67000 Strasbourg, France 

^ Subaru Telescope, National Astronomical Observatory of Japan, 650 North A'ohoku Place, Hilo, HI 96720, USA 

^ Max-Planck-Institut fiir extraterrestrische Physik, D-85740, Garching bei Miinchen, Germany 

* Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 

Received ; accepted 

Abstract. We obtained deep optical imaging of the thermally emitting X-ray bright and radio-quiet isolated 
neutron star RX J1605. 3-1-3249 with the Subaru telescope in 1999 and 2003. Together with archival HST images 
acquired in 2001 these data reveal a proper motion oi n = 144.5 ± 13.2mas/yr. This implies a relatively high 
spatial velocity and indicates that the star is unlikely to be re-heated by accretion of matter from the interstellar 
medium. Assuming that RX J1605. 3-1-3249 is a young (10^-10^ yr) cooling neutron star, its apparent trajectory is 
consistent with a birth in the nearby Sco OB2 OB association at a location close to that derived for RX J1856.5- 
3754 and perhaps also to that of RX J0720. 4-3125. This suggests that the X-ray bright part of ROSAT-discovered 
isolated neutron stars is dominated by the production of the Sco 0B2 complex which is the closest OB association 
and a part of the Gould belt. The B and R magnitudes of the faint optical counterpart did not vary from 1999 
to 2003 at B = 27.22±0.10. Its B-R colour index of -f0.32±0.17 is significantly redder than that of other isolated 
neutron stars and the optical flux lies a factor 11.5 above the extrapolation of the X-ray blackbody-like spectrum. 
The red optical colour reveals the presence of an additional emitting component in the optical regime over the 
main neutron star thermal emission. We also discovered a small elongated Ha nebula approximately centered on 
the neutron star and aligned with the direction of motion. The width of the nebula is unresolved and smaller 
than ~ 0.4" for a length of about 1". The shape of the Balmer emitting nebula around RX J1605. 3-1-3249 is very 
different from those seen close to other neutron stars and should be confirmed by follow-up observations. We 
shortly discuss the possible mechanisms which could give rise to such a geometry. 



1. Introduction 

A decade ago, it was expected that a large number of 
isolated neutron stars (INSs) accreting from the inter- 
stellar medium would populate the ROSAT all-sky sur- 
vey. It is now clear that this numerous population is not 
detected, although a handful of X-ray bright INSs were 
found. These INSs are radio-quiet, display thermal X-ray 
spectra with kT ^ 40-100 eV, undergo little interstellar 
absorption indicating relatively small distances (few hun- 
dred parsecs) and are not associated with any supernova 
remnant (see reviews in Treves et al. 120001 Motch lTOOll and 
Haberl E003|l . Four of them exhibit pulsations with periods 
around 10 s (Haberl |2004 ) . Most of these INSs now have 
identified faint optical counterparts or good candidates 
with B magnitudes in the range of 25.8 to 28.6. Their prox- 
imity and the apparent absence of strong non-thermal ac- 
tivity makes them unique laboratories for testing radiative 
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properties of neutron star surfaces, high gravity and high 
magnetic field physics. This group of INSs could eventu- 
ally bring important constraints on the debated equation 
of state of matter in neutron star interior. 

Grating X-ray spectra of the two brightest ROSAT- 
discovered INSs RX J0720.4-3125 (Paerels et al. .2001|i) and 
RX J1856.5-3754 (Burwitz et aL EUHTl Drake et al. lMi^ 
failed to reveal any spectral feature that could be used to 
derive a gravitational redshift, for instance. Their X-ray 
spectra are usually best fitted by a blackbody-like model. 
Further observations with the EPIC and RGS instruments 
on board XMM-Newton have revealed in several cases the 
presence of broad absorption lines with central energies 
in the 300-500 eV range on top of the blackbody energy 
distribution. These structures have been tentatively inter- 
preted as proton cyclotron lines in magnetic fields of a 
few 10^3 Q (Haberl et al. 1^001 Haberl et aL l^UUl and van 
Kerkwijk et al. l2004|l . A puzzling problem is that none of 
the models fitting the X-ray data can adequately represent 
the optical-UV emission. Hydrogen atmosphere models 
usually overpredict optical flux while the Rayleigh-Jeans 
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tails of the X-ray blackbody model typically fall short by 
a factor « 10 at A « 550 nm (see Pavlov et al. 120021 for a 
recent review). Understanding neutron star surface emis- 
sion properties is therefore a prerequisite to derive reliable 
stellar radii. 

When observed, the X-ray pulsations can be inter- 
preted as due to a non uniform temperature distribution 
over the surface of a young cooling neutron star in strong 
magnetic field conditions or to variable absorption lines. 
The absence of detected pulsations in cases where evidence 
for strong magnetic fields exists is also not well under- 
stood. 

The evolutionary status of these objects, old accret- 
ing neutron stars or relatively young cooling stars can be 
established by various means. A braking or an acceler- 
ation of the rotation period can in principle reveal the 
powering mechanism (Zane et al. I2002|l . This exercise is 
however difficult and limited to pulsating stars. A poten- 
tially more efficient way is to search for proper motion 
using either deep optical imaging or by taking advan- 
tage of the high spatial resolution affordable by Chandra. 
Because of the strong dependency of Bondi-Hoyle ac- 
cretion efficiency on velocity relative to the interstellar 
medium, a high proper motion will naturally exclude the 
accretion scenario, unless high interstellar densities are lo- 
cally encountered. Significant spatial motions have been 
already detected from RX J1856. 5-3754 (Neuhauser"OTjr, 
Walter EOHTl Kaplan et al. .2002,1 and from RXJ0720.4- 
3125 (Motch et al. IMTSjl which in both cases appear to 
exclude the accretion scenario. 

We report here on deep optical imaging obtained at the 
Subaru telescope and Hubble Space Telescope of the op- 
tical counterpart of RX J1605.3+3249 (Motch et al lTU^ 
which is the third X-ray brightest of its class. An optical 
counterpart has been identified by Kaplan et al. I)2()()3a|l . 
Our observations covering a four year interval of time re- 
veal the high proper motion of the object. A small but 
resolved Ha nebula is also detected. 

2. Optical observations 

Optical imaging has been obtained on two occasions with 
the Subaru telescope in 1999 and 2003, and in 2001 by the 
HST. 

Data from the first epoch have been retrieved from the 
SMOKA database at the end of January 2003. These first 
observations were obtained between April 22 and May 15 
1999 using the SuprimeCam instrument (Miyazaki et al. 
I1998|l mounted at the Cassegrain focus. This configura- 
tion which used a 3x2 CCD setup was in operation until 
July 1999 to perform an initial detailed testing of both 
the telescope and instrument. Each ST-002A SITe CCD 
had 2048 x 4096 square pixels with a size of 15/i corre- 
sponding to 0.03" on the sky at the Cassegrain focus. For 
our purpose, we only used CCD # 5 which fully contains 
the field of RXJ1605. 3+3249. Of the 3 filters used, B, 
R and broad Hq , there exist only flat-fields for the R 
band. The Ha filter NA_L651 has a central wavelength of 



652lA and a full width of 325A. Due to relatively large 
central wavelength variations with off axis angle (up to 
100 A), this filter was later decommissioned. Laboratory 
measurements show that the Ha line was in all cases well 
inside the filter band pass but due to the spatial inhomo- 
geneities the exact transmission at Ha remains uncertain. 
Individual exposure times were 1200 s for the Ha and B 
filters and 900 s for the R band. The telescope was moved 
by about 12" in both directions between two consecutive 
exposures. All individual images were corrected for bias 
using over-scan areas and the R band images were further- 
more flat-fielded using well illuminated dome flats. Since 
only one bias exposure exists for these observations, we 
prefered to correct images using mean over-scan values in 
order to not increase photometric and positional errors. 
Individual images were then moved to a common frame 
using a well exposed single reference star located close to 
the neutron star position in order to minimize the effects 
due to the lack of correction for geometrical distortion. 
The images were then stacked using a statistical cosmic- 
ray event rejection criterion and flnally rebinned with a 
0.15" pixel size. 

The HST archive provided second epoch observations. 
We used the unfiltered CCD (50CCD aperture) image ob- 
tained on July 21 2001 which offers the highest sensitiv- 
ity. Individual HST images were geometrically corrected, 
drizzled onto a reference image with a pixel scale half 
of the original one, filtered for cosmic-ray impacts and 
stacked together following the method outlined in Fruchter 
& Hook (2002). These data were presented in Kaplan et 
al. H2003a|l . 

Third epoch images were all obtained on June 8 2003 
with the FOCAS instrument (Kashikawa et al. 12002(1 op- 
erated at the Cassegrain focus. The target was located 
on CCD # 2 which has the best cosmetic quality. We 
obtained series of B and R band exposures with integra- 
tion times of 1200 s and 600 s respectively. Raw images 
were corrected for bias using average offset exposures and 
flat-fielded with dome flats. A dedicated IDL procedure 
then corrected individual images for geometrical distor- 
tion which otherwise could lead to errors of up to 2" at 
the edge of the held of view. All images were moved to 
a common frame using a set of 9 reference stellar-like ob- 
jects and then stacked together using a statistical criterion 
for rejecting cosmic-ray impacts. 

Ground based observations were in general obtained 
in photometric conditions although some thin cirrus may 
have been present in 2003. The Subaru telescope is 
equipped with an atmospheric refraction correction prism 
at the telescope side. The difference in colours between 
the relatively red astrometric reference stars and the blue 
neutron star thus does not yield any signiflcant difference 
in position. 

The log of observations is listed in Table HI FWHM 
seeings are given as measured on the summed images. 



C. Motch et al.; The proper motion of RX J1605. 3+3249 



3 



Table 1. Optical data 
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Table 2. B and R magnitudes of RX J1605.3+3249 

Epoch B R 

April-May 1999 27.25 ± 0.19 26.90 ± 0.15 
June 2003 27.22 ± 0.10 26.90 ± 0.14 



3. Photometry 

During the second Subaru observations, B and R pho- 
tometric standard stars were repeatedly observed in the 
field of SAllO-504 fLandolt IT??^ . However, zero point 
variability by about 0.1 mag rms reveals the probable 
presence of thin cirrus during that night. We thus used 
our photometric CFHT images obtained in 1998 and de- 
scribed in Motch et al. H1999|l to calibrate several B and 
R local photometric reference stars in addition to those 
already mentioned in that paper. Comparing CFHT and 
Subaru 2003 zero points shows that the 2003 images were 
indeed absorbed by ^ 0.15 mag in B and R. Magnitudes 
were computed using a 2-d Gaussian fitting process tak- 
ing into account possible sky background spatial variabil- 
ity. The magnitudes of RX J1605. 3+3249, hsted in Table 
121 are generally consistent with those reported by Kaplan 
et al. l|2UU3all from HST/STIS photometry. Errors quoted 
are one sigma and take into account uncertainties in colour 
transformation and for the 1999 B value errors due to the 
absence of flat-field correction. There is no evidence for 
flux variability on time scales of years. We show in Fig. 
n] the summed B image obtained in 2003 at the Subaru 
telescope with FOCAS. 

Using the best determined value in 2003 we find B— R 
= +0.32 ± 0.17. This colour index is significantly redder 
than that of RXJ1856. 5-3754 (B-R = -0.61 ± 0.13, van 
Kerkwijk & Kulkarni 2001.3) and maybe also RX J0720.4- 
3125 (B-R = -0.3 ± 0.4, Kulkarni & van Kerkwiik lTM5|l . 

4. Proper motion study 

Comparing positions obtained by three different instru- 
ments with different band passes requires some caution. 
We estimated the errors resulting from uncorrected geo- 
metrical distortions by measuring the dispersion in posi- 
tion of a number of bright stellar-like objects on individual 
Subaru images. For the 1999 observations, the rms error is 
12 mas in the B filter (4 objects) and 7 mas in the Ha and 
R filters (5 objects). Geometrical correction of the 2003 




Fig. 1. The summed B band image obtained in 2003 with 
FOCAS 



FOCAS frames brings down the dispersion from about 12 
mas to 4.6 mas. Therefore, the absence of available distor- 
tion mapping for the Suprime-Cam instrument mounted 
at the Cassegrain focus has little impact. However, the 
geometrical correction of FOCAS images indeed improved 
the astrometric quality. 

We searched for proper motion using the sum of the 
B and R band images in 1999, the HST 50CCD image 
and the B band FOCAS image in 2003 which provide the 
best compromise between sensitivity and image quality. A 
group of astrometric reference objects were selected on the 
HST image for their stellar-like profile or very small spatial 
extension. We used a 2-d Gaussian fitting process allowing 
for sky gradient to measure accurate positions. Centering 
accuracies are listed in Table 01 The Subaru 1999 and 
2003 images were scaled and rotated to the HST frame 
using 9 and 7 such astrometric reference stars respectively. 
After transformation to the HST frame, reference objects 
had mean position residuals of the order of 16 mas and 
25 mas for the Surprime-Cam and FOCAS, respectively. 
These residuals are larger than expected from centering 
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Fig. 2. Positions of RX J1605. 3+3249 for the 3 epochs 
(fiUed squares) plotted in the reference pixel frame of the 
HST/STIS image. The best fit proper motion solution 
(/i = 144.5 ± 13.2mas/yr) is shown as open squares 

position errors only and probably reflect some remaining 
low level large scale field distortion. The Subaru positional 
accuracies transformed to the HST reference frame take 
these systematic errors into account. 

A proper motion solution expressed in terms of to- 
tal motion and direction was then fitted to the three 
recorded positions of RX J1605. 3+3249 using a mini- 
mum algorithm. The best fit shown in Fig. |21 is 
obtained for /i = 144.5 ± 13.2mas/yr with a direc- 
tion of motion of —9.86° ± 5.65° eastward from north. 
This corresponds to cos(5)/Xq,=— 24.7±16.3mas/yr and 
/X5=142.4±15.4mas/yr. All errors are quoted here for a 
90% confidence level on one parameter. The detected 
proper motion is thus highly significant and overall mea- 
surement and systematic errors appear negligible (see Fig. 
O. The part of the apparent displacement due to solar 
motion should be small, of the order of 5 mas in right as- 
cension and 11 mas in declination at 100 pc. 

Most of the proper motion of RX J1605.3+3249 is in 
galactic longitude (/i; — 137 mas/yr). The source is how- 
ever slowly moving away from the galactic plane with fif, 
= 37 mas/yr. The unknown distance and radial velocity 
of the neutron star does not allow us to derive the pos- 
sible position of its birth place with very good accuracy. 
In order to constrain the origin of RX J1605. 3+3249 we 
assumed a maximum age of ^ 10^ yr, compatible with 
possible cooling times of a young neutron star with kT 
~ 100 eV (see below, Yakovlev & FethickEOMIl, distances 
in the range of 100 pc to 500 pc and radial velocities in 
the range of ± 700kms~^. Considering the short flight 
time, we also neglect the effects of the galactic potential. 
The only nearby OB association overlapping with some 



motion of RX J1605. 3+3249 




Fig. 3. Contours of the Subaru observations plotted over 
the HST 50CCD image 



Table 3. Centering accuracies (mas) 

Sub. 1999 HST 2001 Sub. 2003 
Reference stars 8 3 9 
RXJ1605. 3+3249 25 4 10 



backward trajectories is the upper Scorpius-Centaurus 
region belonging to the Sco 0B2 complex. Amazingly, 
Walter & Lattimer H2UU2|I derive a similar birth place for 
RX J1856. 5-3754 in the upper Scorpius while Motch et al. 
(2003) speculate on a possible birth place in the lower 
Centaurus-Crux part of the same Sco 0B2 complex for 
RXJ0720.4-3125. 

5. An Ha nebula 

We show in Fig.0|the Ha image obtained in 1999. An ex- 
tended and elongated source is clearly seen at the position 
of RXJ1605. 3+3249. The emission is detected on individ- 
ual images whose centers were shifted by several hundreds 
of pixels and is thus unlikely to be of instrumental ori- 
gin or caused by the absence of flat-field correction. The 
direction of largest extent of the nebula is close to that 
of the apparent trajectory. Ha nebulae have been discov- 
ered around a small number of isolated neutron stars (see 
Chatterjee & Cordes .2002. for a recent review). They have 
in general a well defined arc-like shape extending several 
tens of arcseconds on the sky which can be interpreted as 
due to a bow shock created by the interaction of the high 
velocity pulsar wind with the interstellar medium. In the 
case of RX J1856. 5-3754 Ha emission could also be due to 
X-ray ionization (van Kerkwijk & Kulkarni ■2001b)l . The 
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bright cylindrical head of the Guitar nebula produced by 
the radio pulsar B2224+65 (Cordes et al.UM^. although 
significantly larger, could be a scaled-up analogue of the 
nebula seen here. However, we shall see below that this is 
in fact unlikely and that the shape and probably also the 
origin of the nebula detected around RX J1605. 3+3249 are 
unique. 

We tried to constrain the actual size and orientation of 
the Ha emission by fitting to the observed image a model 
consisting of a cylinder convolved with the point spread 
function. Owing to the poor spatial resolution available, 
we assumed a simple uniformly emitting bar shape charac- 
terized by its length, width and orientation. The Gaussian 
parameters of the point spread function were derived from 
the nearby star "C" which is unresolved in the HST im- 
age and is a likely dwarf M star (Motch et al. 11999(1 . We 
used a 2-d minimum fitting process taking into account 
CCD readout noise and counting statistics. The nebula 
seems slightly wider than would be expected from the PSF 
alone. However, the resulting width of the bar depends on 
whether the background is locally estimated or computed 
from a wider area as well as on the overall size of the fitted 
image. The possible width ranges from 0", i.e. unresolved, 
to « 0.4". The best fit length is « 1.1" with an inclina- 
tion angle of — 4.5°± 5.5° (Ict error) eastward from north, 
consistent with the apparent direction of motion. 

The Ha nebula is approximately centered on the neu- 
tron star (see Fig.EJ with a slightly larger extent (« 0.6") 
toward the direction of motion. A total of 7 stellar like ob- 
jects were used to move the Ha image onto the combined 
R+B image with a rms positional accuracy of 12 mas. The 
off center position of the neutron star with respect to its 
associated nebula is thus significant. 

In the absence of flat-field and filter calibration at the 
position of RX J1605. 3+3249 we can only very roughly es- 
timate the Ha flux from the nebula. We assume a mean 
filter transmission of 80% similar to that of the Ha filter 
currently mounted on FOCAS and a total calibration er- 
ror of 20%. The integrated flux from the nebula amounts 
to 2.4+0.6 lO"'^ photon s'^ cm"^ or 7.2+1.8 10"^^ erg 
s~^ cm~^. Although much less extended than that of 
RXJ1856. 5-3754 (van Kerkwijk & Kulkarni l20(nb|l . the 
Ha surface brightness of RXJ1605. 3+3249 seems about 
10 times brighter (for a 0.4" width) allowing its detection 
in a Ih exposure only on an 8-m class telescope. 



6. Discussion 

6.1. Proper motion 

RX J1605. 3+3249 displays the second largest proper mo- 
tion of the three objects of this class for which measure- 
ments exist; RX J1856. 5-3754 has ^ = 333 ± Imas/yr 
(Neuhauser 2001 , Walter IMJTl Kaplan et al. I^02j) and 
RXJ0720.4-3125 has n = 97±12 mas/yr (Motch et al. 
nnn^ . For RXJI6O5.3+3249, a proper motion of = 
144.5 ± 13.2 mas/yr corresponds to a transverse velocity 




Fig. 4. The Ha image obtained in 1999 with the Suprime- 
Cam camera. The arrow marks the neutron star position 




Fig. 5. Overlay of the Ha image on the R+B image ob- 
tained in 1999 with the Suprime-Cam camera. One pixel 
corresponds to 0.15" on the sky 



of Vt ~ 68 X (d/lOOpc) kms ^ with a small component 
due to solar motion of ~ 5.5x (d/100pc)~^ kms^^. 

A bolometric luminosity of 1.1 x 10"^^ (d/lOOpc)^ erg 
s^^ (Motch et al. I1999|l derived from the blackbody fit 
to the ROSAT data, if powered by accretion from the in- 
terstellar medium, would require a local particle density 
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^ 170 (d/lOOpc)'^ cm^'^. This furthermore assumes that 
there is no radial velocity component. Such ISM densi- 
ties are much larger than the average density toward the 
source; n ~ 0.3 (d/100pc)~^ cm""^, estimated from the 
hydrogen column density required by the blackbody fit of 
the ROSAT X-ray spectrum (Nh'-^ 10^° cm'^, Motch ct al. 
I1999p . Following the arguments in Motch et al. (|2003,) we 
can conclude that like RX J0720.4-3125, RX J1605.3-h3249 
is unlikely to be powered by accretion from the ISM un- 
less extremely rare conditions are met. Similar conclusions 
were reached by van Kerkwijk & Kulkarni (2GGlb) for 
RXJ1856. 5-3754 based on the modeling of its associated 
Ha nebula. 



6.2. Is Sco 0B2 the main provider of the brightest 
radio quiet INSs 7 

RXJ1605. 3-1-3249 is thus the third case of X-ray bright 
and radio quiet INS for which the accretion scenario can 
be excluded with relatively high confidence. This confirms 
the conclusion that the thermally emitting neutron stars 
detected by ROSAT are in general not powered by accre- 
tion from the interstellar medium but are rather young 
cooling neutron stars with typical ages ranging from 10^ 
to 10^ yr. The absence of accurate distance estimates, ra- 
dial velocity information and the still relatively large er- 
rors on proper motion direction of RXJ1605. 3-1-3249 and 
RX J0720. 4-3125 do not allow us to draw unambiguous 
conclusions on the birth place of these neutron stars. 
There are however relatively strong hints that the three X- 
ray brightest ROSAT discovered INSs, RX J1856.5-3754, 
RXJ1605. 3+3249 and possibly RX J072G.4-3125 were all 
born in the Sco 0B2 association. With a mean distance of 
120-140 pc (de Zeeuw et al. ll999|l this rich OB association 
is the closest to earth and belongs to the large structure of 
young stars known as the Gould belt. Since most OB stars 
are found in associations, we expect that statistically at 
least, the fact that a neutron star trajectory crosses an OB 
association at a time roughly consistent with the cooling 
time is evidence for its birth place. Some other middle- 
aged (age < 4.25 Myr) nearby (d < 350 pc) neutron stars 
have tentative birth places, Geminga possibly in Orion 
OBI (Caraveo et al.!TW5|, PSR J1932+1059 also in Sco 
0B2 (Hoogerwerf et al.^EI- The revised distance of PSR 
B0656+14 (d - 290 pc Brisken et al. lMI^ could also sug- 
gest a birth place in Orion OBI. The tangential velocities 
of the X-ray bright and radio-quiet INSs have relatively 
modest values for neutron star standards, 220±60 kms~^ 
for RXJ1856. 5-3754 (Kaplan et aLHHH, 50 (d/lOOpc) 
kms-i for RXJ0720.4-3125 (Motch et al. EPj and 68 
(d/lOOpc) kms^i for RX J1605.3+3249. Their velocities 
are within the range observed for radio pulsars and, as- 
suming distances of 200 pc for RX J0720.4-3125 (Motch 
et al. lMTSjl and ^ 300 pc for RX J1605.3+3249 (Kaplan 
et al. I2()03a|l . suggest that they belong rather to the slow 
part of the pulsar population (Arzoumanian et al. I2()()2fl . 



Although still based on small number statistics, it 
seems that our current proximity to Sco 0B2 may be the 
reason why we can detect in relatively large numbers X-ray 
bright and radio quiet neutron stars at small distances. In 
contrast, Orion OBI, the other major OB association in 
our neighborhood may have provided Geminga and PSR 
B0656+14 but so far no ROSAT INS. Matching the flight 
time from Orion (d ^ 500 pc) and the cooling time im- 
plies large radial velocities (~ 700 kms~^). Such large 
radial velocities compared to tangential (120 kms~^ for 
Geminga) could be understood as due to a selection effect 
since we only detect in X-rays the neutron stars born in 
Orion which are young, hot and fast enough. The possibil- 
ity is thus open that radio quiet and hot neutron stars may 
preferentially be born with low kick velocities. However, 
before such conclusions can be drawn, proper motions and 
distances of fainter objects must be measured. 

6.3. Optical to X-ray energy distribution 

We show in Fig.lHlthe overall X-ray to optical distribution 
of RX J1605.3+3249. The model describing the X-ray con- 
tinuum is a simple blackbody fit to EPIC pn data obtained 
in small window mode through the thick filter during rev- 
olution 589. We preferentially used these data because of 
their low pile-up fraction and high quality of the calibra- 
tions available for this instrumental setting compared to 
others, especially the thin filter which suffers from some 
transmission uncertainties at low energies. The blackbody 
fit gives a of 196 for 174 dof. Adding a Gaussian line 
with a fixed energy of 0.45 keV (van Kerkwijk et al. l2()()4|l 
significantly improves the fit, reducing the cit 167.5 for 
172 dof. The best fit Gaussian equivalent width is 29.5 eV 
with a width a ^ 180 eV. For simplicity we plot in Fig. |B1 
the best blackbody-only fit with kT = 98.9 eV and R°° 
0.9 (d/100pc)km. 

As for all X-ray dim optically identified INSs, the 
optical continuum of RX J1605. 3+3249 lies above the 
Rayleigh-Jeans tail of the blackbody-like spectrum re- 
sponsible for the X-ray energy distribution. This has been 
sometimes used as an argument in favour of the exis- 
tence of a second cooler thermal component which would 
hardly contribute in X-rays but could account for most 
of the optical emission (see Pavlov et al. I2002|l . In the 
case of RX J1856. 5-3754 , for which there exists an ac- 
curate distance estimate based on a measured parallax, 
the low temperature implies a larger radius of the emit- 
ting area more consistent with equations of state with 
conventional hadronic matter (Burwitz et al. I2003|l . The 
two-blackbody model can naturally represent the strong 
emissivity anisotropy induced by a large dipolar magnetic 
field. In this framework, the hot thermal component would 
represent radiation emitted from the polar regions while 
the cooler one would be emitted by the equatorial regions. 
However, in the case of RXJ1605. 3+3249 (van Kerkwijk 
et al. Enni) and even more acutely for RX J1856.5-3754 
(Ransom et al. I20()2|l , the absence of pulsations down to 
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low pulse fractions puts severe constraints on the geomet- 
rical configuration and amplitude of gravitational deflec- 
tion. Standard hydrogen atmosphere models, while usu- 
ally fitting correctly the X-ray energy distribution, have 
the common problem to largely overpredict the optical 
continuum. Thin hydrogen atmospheres such as proposed 
by Motch et al. iPMjl for RX J0720.4-3125 can account 
for the optical to X-ray energy distribution with a single 
component but have the drawback of assuming the pres- 
ence of an ad hoc amount of light elements on the surface 
of the neutron star. A high surface reflectivity in the X-ray 
energy range such as in the case of condensed atmospheres 
might also solve the X-ray - optical discrepancy (Burwitz 
et al. lMiSji . 

The B-R colour index of RX J1605.3+3249 (B-R 
= -1-0.32 ± 0.17) is significantly redder than that of 
RXJ1856. 5-3754 which essentially follows a F(A) ~ A"^ 
continuum. Our B band flux is consistent with that mea- 
sured in the wide HST/STIS 50CCD band pass whereas 
our R band flux is somewhat above that computed for 
the HST/STIS F28X50LP long pass filter. At present, the 
reasons for this small discrepancy (formally a 2.1 cr effect) 
are unclear. Due to flux calibration uncertainties, the con- 
tribution of the nebula detected in the broad Ha filter to 
the R band flux is difficult to estimate. The fact that the 
nebula remains undetected in the R band images suggests 
that the effect is not strong. The very small Nh required 
to fit the X-ray data (Nr < 2.2 10^" cm-^, van Kerkwijk 
et al. 12004(1 . corresponds to a maximum E(B-R) of about 
0.09, too small to account for the red colour index if an 
intrinsic Rayleigh-Jeans spectrum is assumed. 

The improved spectral resolution clearly shows that 
the energy distribution in the optical range does not fol- 
low a Rayleigh-Jeans tail and reveals the presence of 
an additional component. A similar situation exists in 
RX J0720. 4-3125 where a power law with an index slightly 
softer than Rayleigh-Jeans is required to fit the opti- 
cal/UV continuum (Kaplan et al. I2003b)l while still keep- 
ing relatively blue overall colours. Generally, the B— R 
colour index of RXJ1605. 3-1-3249 appears similar to that 
of two middle aged neutron stars, Geminga (B— R ~ 
0.2, Bignami et al. 1993j and PSR B0656-f-14 (B-R 
0.62±0.11, Komarova et al. I2003II . In these two pulsars, 
however, magnetospheric activity exists in the X-ray en- 
ergy range (Zavlin & Pavlov 2004). Mignami et al. (1998) 
find evidence for a peak at ~ 6OOOA which they interpret 
as the possible signature of cyclotron emission. For PSR 
B0656-f 14, a single power law component can in principle 
explain both the hard X-ray tail and the excess of op- 
tical and near-infrared emission above a Rayleigh-Jeans 
component (Pavlov et al. 2002). Assuming that a similar 
mechanism is at work in RXJ1605. 3-1-3249, we computed 
that a single powerlaw component with a photon index of 
about 1.5, similar to that found in PSR B0656-I-14, could 
account for both the 2-10 kcV EPIC pn flux 1 a upper 
limit (1.8 10~^^ergs~"'^ cm~^) and for the B to R band 
flux. The magnetospheric activity of RXJ1605. 3-1-3249 
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Fig. 6. The spectral energy distribution of 
RX J1605.3-f3249. Dotted hue: unabsorbed blackbody 
fitting the X-ray XMM-Newton EPIC pn data (see text). 
Thick line: corrected for absorption. The two data points 
represent our B and R band measurements. The dashed 
curve is a power-law spectrum of photon index 1.53 and 
the dashed-dotted curve gives the sum of the blackbody 
and powerlaw spectra. Zoomed frame: optical energy 
distribution based on HST and Subaru measurements 



may thus well be a scaled down version of that of PSR 
B0656-M4. 

6.4. Nature of the Ha nebula 

Ha nebulae have been discovered around about six INSs, 
ordinary or millisecond pulsars. In all cases an arc-shaped 
emission extending over several tens of arcseconds is vis- 
ible with the main nebula axis roughly aligned with the 
direction of proper motion (Chatterjee & Cordes EHEI) • 
Supersonic pulsar motion can create a bow shock neb- 
ula when the high velocity wind which carries away most 
of the rotational kinetic energy interacts with interstellar 
medium. The shock surface is located where ram pres- 
sure balance between the pulsar wind and the interstel- 
lar medium takes place. Ha emission then mainly arises 
from coUisional excitation of neutrals by shocked electrons 
and protons. In many cases, the shape of the nebula can 
be reasonably well fitted by the model of Wilkin H199fi|l 
which assumes an isotropic wind and rapid radiative cool- 
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ing. Some asymmetric nebulae provide evidence toward 
density gradients in the ISM, pulsar wind anisotropy or 
both (see e.g. Gaensler et al. 12002). 

No such arc-shaped nebula is seen around 
RXJ1605. 3+3249. Our relatively low detection sen- 
sitivity could explain part of the observed differences. 
However, in the known bow shock nebulae Balmer emis- 
sion usually shows the highest surface brightness ahead of 
motion. At the apex the nebula extends perpendicular to 
the trajectory thus tracing the shape of the front-shocked 
region. In the present case, the Ha emission appears 
elongated in the direction of motion. Such a geometry 
cannot be matched by any isotropic wind model since 
one would naively expect a larger extent in the direction 
perpendicular to the motion than in the direction of 
motion. In the absence of detailed modeling, it is unclear 
whether pulsar wind beaming could account for the 
observed geometry. 

If no strong wind blows from RX J1605. 3+3249 Ha 
emission could still arise from an X-ray ionized nebula. 
Blaes et al. (1995) have shown that ionization of the local 
interstellar medium by the X-ray and UV radiation field of 
the hot neutron star combined with proper motion would 
produce cometary HII regions. Their properties depend 
on the details of the ionization, recombination, heating 
and cooling rates (see van Kerkwijk & Kulkarni I2001b| 
for a detailed modeling including dynamical effects). The 
isotropic ionizing field should create a spherical head neb- 
ula centered on the neutron star with an extent in the 
direction perpendicular to motion at least as large as in 
the direction of motion, similar to the case of bow shock 
nebulae. Also, because of the rather long cooling time, 
most of the Ha emission will occur in the wake of the neu- 
tron star, again at variance with the observation. Scaling 
the observed Ha surface brightness to that of RX J1856.5- 
3754 (van Kerkwijk & Kulkarni 1200 lb|l would imply local 
densities of the order of nn ~ 9 (d/100pc)~^ cm~^, much 
larger than expected for normal ISM densities. Such high 
densities would however yield a cooling time of TcooI 
6(d/100pc) yr roughly comparable to the nebula cross- 
ing time Tcross ~ 4yr as estimated from its angular size 

600 mas ahead) and observed proper motion. 

The geometry of the nebula is in principle consistent 
with that expected from a bipolar jet. Such jets have been 
observed in X-rays from very young neutron stars, the best 
documented cases being those of the Crab (Weisskopf et 
al. rWK^ and of the Vela pulsar (Helfand et al. IMIUl . 
Their power-law X-ray spectra are interpreted as due to 
optically thin synchrotron emission from ultra-relativistic 
particles. Obviously, some other mechanism must be in- 
voked to explain Ha emission at redshift < 0.03 (the max- 
imum allowed by the filter band pass) from such a jet. 
One possibility could be the interaction of the beamed 
particle flow with the interstellar medium. If the main 
direction of the Ha nebula would indeed trace the spin 
axis, RX J1605. 3+3249 would be another case of appar- 
ent alignment of spin axis and proper motion direction 
(Lai et al. i2001)l . 



Alternatively, if the main axis of the nebula would ex- 
actly follow the trajectory, as would be the case for an 
X-ray ionized or bow shocked nebula unperturbed by the 
interstellar medium, the apparent proper motion direction 
should not be exactly aligned with the nebula as a result of 
solar motion. In principle, this effect could be used to de- 
rive independent estimates of the distance of the neutron 
star. However, the accuracy with which we can determine 
the angle difference 9 = 4.6°± 6.3° is not yet large enough 
to be useful. Nevertheless, the sign of the angle difference 
between the nebula and the apparent proper motion is 
that expected due to the different positions of the Sun 
between 1999 and 2003. Using the Sun velocity param- 
eters determined by Hipparcos (Dehnen & Binnev [l998|l 
we compute an angle difference 9 = 1.3°(d/100pc)~^ . 

Considering the absence of calibration of the instru- 
mental Ha setting, it is clear that the astonishing geom- 
etry of the nebula needs to be confirmed and specified 
by other observations before extensive conclusions on its 
nature can be drawn. 

7. Conclusions 

Our deep Subaru imaging reveals the high proper mo- 
tion of RX J1605. 3+3249 and thus confirms the optical 
identification proposed by Kaplan et al. H2003a|l . The cor- 
responding large space motion most probably excludes 
accretion from the interstellar medium as the source of 
the radiated X-ray energy. Similarly to RX J1856. 5-3754 
and RX J0720.4-3125, RX J1605. 3+3249 is likely a cooling 
neutron star with an age younger than ^ 10® yr. Evidence 
for residual magnetospheric activity comes from the rel- 
atively red B— R index which indicates the presence of 
an excess of red light above the expected Rayleigh-Jeans 
tail, a picture comparable to that seen in Geminga or in 
PSR B0656+14. A small size Ha nebula is detected ap- 
proximately centered on the neutron star and extended 
along the direction of the trajectory. The lack of photo- 
metric calibration and flat-fielding does not make it possi- 
ble to draw stringent conclusions on its nature. However, 
the small size and geometry of the nebula are unique and 
no standard bow-shock or X-ray ionisation model can ac- 
count for these peculiarities in a straightforward manner. 
Forthcoming deeper observations should allow us to con- 
firm the unusual features detected in 1999. Finally, we find 
that RXJ1605. 3+3249 may have been born in the same 
nearby Sco 0B2 OB association as RXJ1856. 5-3754 and 
also perhaps RX J0720. 4-3125. This could imply that the 
bright part of the observed LogN-LogS curve for X-ray dim 
isolated neutron stars is dominated by the production of 
this single OB association. 
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